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We demonstrate enhancement of the dispersive frequency shift in a coplanar waveguide resonator 
induced by a capacitively-coupled superconducting flux qubit in the straddling regime. The mag- 
nitude of the observed shift, 80 MHz for the qubit-resonator detuning of 5 GHz, is quantitatively 
explained by the generalized Jaynes-Cummings model which takes into account the contribution of 
the qubit higher energy levels. By applying the enhanced dispersive shift to the qubit readout, we 
achieved 90% contrast of the Rabi oscillations which is mainly limited by the energy relaxation of 
the qubit. 
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High-fidelity readout is essential for quantum informa- 
tion processing. In the architectures of circuit quantum 
electrodynamics based on superconducting circuits [![, 
dispersive readout using a state-dependent shift x °f the 
cavity resonance is commonly used @, Q . In this scheme, 
it is important to choose a proper value of the cavity de- 
cay rate k in order to maximize the signal-to-noise ratio. 
Large k is favorable in terms of the information flow rate, 
but unfavorable in terms of the state distinguishability 
As shown in Ref. 4, under a fixed number of photons in 
the resonator, the signal-to-noise ratio is optimized by 
k = 2\ for a given x- Thus, large x is advantageous 
because it allows large optimized k. Achieving large x is 
also important for applications such as conditional phase 
gate Q and photon number resolving [ff]. 

The magnitude of x can be increased by making the 
coupling between the qubit and the cavity larger or de- 
creasing the detuning between the cavity and the qubit, 
but both of them are at the price of reducing the up- 
per limit of the qubit relaxation time through the cavity, 
so-called Purcell effect 0. Koch et al. pointed out an- 
other possibility of utilizing the third level of transmon 
qubits [8J. The key requirement is to realize the condi- 
tion wrji < w r < uj 12, where w r is the cavity resonant fre- 
quency and uiij is the qubit transition frequency between 
i-th and j-th levels. Achieving this condition, which they 
call a 11 straddling regime" , enables the cooperative inter- 
play between 0-1 and 1-2 transitions, giving rise to the 
enhancement of x- However, in order to be in the disper- 
sive regime, it is also required that the coupling strength 
g between the qubit and the cavity is much smaller than 
the qubit anharmonicity, |wi2 — ^oi\, which is only a few- 
hundred MHz for transmons. 

In this Letter, we demonstrate the straddling regime 
in a superconducting flux qubit coupled to a coplanar 



waveguide (CPW) resonator, and obtain 2% as large as 
27rx 80 MHz with a detuning of 5 GHz at the optimal bias 
point of the qubit. The large anharmonicity of the flux 
qubit allows g/27r of a few- hundred MHz while staying 
in the deep dispersive regime. The flux qubit is coupled 
to the resonator by a capacitance instead of an induc- 
tance typically used @. We find that the flux qubits 
can be strongly coupled to the resonator by the capaci- 
tance [10( , and actually the capacitive coupling is impor- 
tant for the observed large x because it produces large 
matrix clement for 1-2 transition. The capacitive cou- 
pling is also advantageous in terms of the readout back- 
action, as demonstrated in a transmon another dc- 
charge- insensitive device. Using the large dispersive shift, 
we achieved 90% contrast of the Rabi oscillations which 
is mainly limited by the energy relaxation of the qubit. 

Figure 1(a) shows the circuit diagram of the device 
we measure. The flux qubit is a conventional three- 
Josephson-junction flux qubit, in which one junction is 
made smaller than the other two by a factor of a. One of 
the two islands defining the smaller junction is coupled 
to the center conductor of the CPW resonator via the 
coupling capacitance C c , while the other is connected to 
the ground. The Hamiltonian of the system consists of 
the qubit part, the resonator part, and the coupling part, 
that is, 

H = H q + H x + H c . (1) 
Each term of the Hamiltonian is given as below [ll[ : 
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FIG. 1. (Color online) Circuit diagram and sample images, 
(a) Circuit diagram of the device, (b) Optical image of 
the device. The meandering trace in the center is a half- 
wavelength CPW resonator made of 50-nm-thick Nb. (c) In- 
terdigitated capacitance between the microwave feedline and 
the resonator, which is designed to be 15 fF. The gap of the 
capacitance is 2-/im wide, (d) Magnified image of the qubit 
part. The qubit made of Al is located in the gap between 
the center conductor and the ground plane of the CPW res- 
onator. One of the islands containing the smaller junction 
is made larger (a white rectangle in the center) so that it 
has large enough capacitance C c to the center conductor of 
the resonator, (e) Scanning electron micrograph of the three- 
junction flux qubit. 



where £j=/ $o/27r, <t> Q =h/2e, / = $ /*"> E c =e 2 /2Cj, 
E T =%IJUC V , V = (l + a)(l+7) + /3, W = a(l + j)+/3, 
X = (1 + 2a)(l + 7) + 2/3, Y = 1 + 2a + 2/3, /3=C C /Cj, 
and 7=C C /C r . Here, <& is an external flux, and I and Cj 
are the critical current and the capacitance of the larger 
Josephson junction of the qubit, respectively Si and Ui 
(i=l, 2) are the phase difference across the larger junc- 
tion, and its conjugate variable representing the charge 
number, respectively. The annihilation (creation) oper- 
ator of a photon in the resonator is denoted as a (<v). 
Finally, L r and C r are the equivalent inductance and ca- 
pacitance of the resonator, respectively. 

Figure 1 (b) shows an optical image of the device. The 
half- wavelength CPW resonator is etched out of a 50-nm- 
thick Nb film sputtered on an oxidized high-resistivity 
silicon wafer. The fundamental-mode frequency ui I /2ir 
is 10.656 GHz at / = 0.5, and the Q factor is 650, 
which is limited by the coupling to the feedline via the 
input capacitance C; n [Fig. 1(c)] [l2j]. The flux qubit 
[Fig. 1(e)] was fabricated by standard shadow evapora- 
tion technique, where two Al layers separated by AI2O3 
are deposited from different angles. The qubit is located 
at one end of the resonator, where the electric field is 
the maximum. The coupling capacitance between the 
resonator and the qubit C c is designed to be 4 fF. The 
qubit-control line couples to the flux qubit via a mutual 




FIG. 2. (Color online) Spectroscopy of the coupled system, 
(a) Phase of P as a function of the flux bias and the probe 
frequency. The vacuum Rabi splittings are clearly observed at 
/ = 0.488 and 0.512. (b) |r| as a function of the flux bias and 
the qubit excitation frequency. Annotations on the energy 
levels denote the corresponding excitation from the ground 
state; |i)q|j)i indicates that the qubit and the resonator are 
in the i-th and j-ih states, respectively. The dashed curves 
represent the fit by Eq. fT]). The same data as in (a) is shown 
in the blue box for comparison. There is no data in the regions 
filled with the meshes. 



inductance of ~0.1 pH. The sample is mounted on a di- 
lution refrigerator and cooled to ~10 mK. 

To confirm the coupling between the qubit and the 
resonator, we first measured the reflection coefficient T 
of the resonator using a vector network analyzer. Fig- 
ure 2(a) shows the phase of T as a function of the flux 
bias and the probe frequency uj p /2tt. The power applied 
to the resonator is —138 dBm, corresponding to 0.02 pho- 
tons in the resonator. Vacuum Rabi splittings are clearly 
observed at the points where the qubit energy Hujqi is 
equal to hu) T , indicating that the strong coupling regime 
is achieved. The observed splitting size is 460 MHz. 

In order to measure the energy levels of the qubit in 
the wider frequency range, we applied another microwave 
at the frequency of u> c /2t: to the control port continu- 
ously. We measure |r| at fixed ui p of 2irx 10.656 GHz, 
while sweeping u c and the flux bias. Changes in |r| is 
observed when the qubit is excited, thus revealing the 
energy band structure. Figure 2(b) shows |r| as a func- 
tion of uj c and the flux bias. As shown in the figure, 
each of the observed energy levels can be assigned to 
a particular excitation of the coupled system from the 
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ground state. The resonances observed at around 18 
and 20 GHz are probably due to resonant modes in the 
sample package. By fitting the data to the Hamilto- 
nian [Eq. (1)], we obtained the following circuit param- 
eters: Ej/h=U8A GHz, E c /h=3.268 GHz, a=0.6106, 
and C c =4.079 fP. All of these values are quite consistent 
with our design. From these parameters, 0-1 and 1-2 
transition frequencies of the qubit at / = 0.5 are calcu- 
lated to be uj 01 /2tt =5.513 GHz and oj 12 /2-k =14.54 GHz, 
respectively. Note that w r is almost right at the middle 
of those frequencies. We come back to this point later. 

Next, we measured the dispersive shift of the cav- 
ity resonance depending on the qubit states by using 
pulsed microwaves and heterodyne detection of the re- 
flected probe signal [Fig. 3(b) inset; see also Ref. 11]. 
The qubit is biased at the optimal point (/ = 0.5) where 
the dephasing due to the flux noise is minimal. When 
the qubit is prepared in the ground state, i.e., for the 
case without a qubit excitation pulse, we observe the 
cavity resonance at 10.656 GHz. It is seen in Fig. 3(a) 
as a dip in the amplitude of the normalized reflection 
coefficient T' = (T — r on )/|r o ff — r on |, where r on and 
r o ff are the reflection coefficients obtained on- and off- 
resonance of the cavity, respectively. For the latter, the 
value of r at uj p /2tt = 10.52 GHz are used as a refer- 
ence. When we apply a 7r-pulse to the qubit, on the 
other hand, a new dip appears at lj p /2it = 10.576 GHz, 
while the one at 10.656 GHz is diminished significantly. 
The magnitude of the observed dispersive shift (2|x|) is 
27rx80 MHz. Since the dispersive shift corresponds to 
the difference in the single-photon excitation energy of 
the resonator depending on the qubit states |0) or |1), 
we can estimate 2\x\ from the energy band calculation 
using the circuit parameters obtained from the fitting. 
The calculated 2\\\ is 27rx71.5 MHz, which is close to 
our observation. Presently, we do not fully understand 
the reason for the ~10% discrepancy, though we sus- 
pect that the nonlinearity of the resonator exerts some 
influence on the measured dispersive shift [l2|]. Note, 
however, that both the observed and the calculated 2|x| 
are much larger than that expected from the standard 
Jaynes-Cummings Hamiltonian, in which the qubit is as- 
sumed to be an ideal two-level system and x = 2 /A, 
where g and A are the coupling strength and the detun- 
ing between the qubit and the resonator, respectively. In 
the present case, 2g is 27rx460 MHz [Fig. 2(a)] and A 
is 2vrx 5.143 GHz, which leads to 2|x| of 2vrx20.6 MHz. 
This value becomes even smaller if we take into ac- 
count a contribution of the Bloch-Siegert shift, namely 
2\ X \ = 2g 2 /A-2g 2 /(uj I + ujoi) = 2ttx14 MHz [lj. These 
facts indicate non- negligible contribution from higher en- 
ergy levels of the flux qubit. Below we elaborate this 
point. 

Following the work on transmon qubits [8j, we con- 
sider the generalized Jaynes-Cummings model and take 
into account the effect of the higher energy levels. The 
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FIG. 3. (Color online) Dispersive shift and Rabi oscillations 
observed at /=0.5. (a) Amplitude of the normalized reflec- 
tion coefficient F' as a function of lj p . The blue curve is taken 
when the qubit is in |0) state, and the red curve is when the 
qubit is in |1) state which is prepared by a 7r-pulse of 3-ns 
long. The shift of the dip frequency represents the dispersive 
shift. The residual dip at cj p /27t=10.656 GHz for the qubit 
in |1) state can be explained by the energy relaxation of the 
qubit. (b) High contrast Rabi oscillations. |r'| is plotted 
as a function of the control pulse length. The blue dots are 
the data and the red curve is a fit by a sinusoidal function 
with exponential decay. The inset depicts the applied pulse 
sequence. After the control pulse with the length of At is ap- 
plied, we wait for td before acquiring the readout signal. The 
length of the readout pulse was 400 ns, and the IF frequency 
was 50 MHz. The signal was sampled at 1 GS/s, and a 100- 
ns-long time trace was used to extract the amplitude and the 
phase. 

effective dispersive shift in the present system is given by 
1 N 

XeS = X01 - XlQ + 2 ^(Xjl ~ Xl 3 ~ X 3° + ( 5 ) 
i=2 

where Xij = \9ij\ 2 /(^ij -w r ), gij = (i\H c \j), Wy=w,- -u)i, 
and N is the number of qubit states taken into account. 
Here, %uii represents the energy of the i-th qubit state. 
Unlike the case for a transmon, matrix elements gty can 
be large not only for neighboring levels, i.e., i = j±l, but 
also for other pairs of levels because of large anharmonic- 
ity of the flux qubit. It is also worth mentioning that 
in the above formula, the contributions of the counter- 
rotating terms (xij with i > j) are taken into account. 
As shown below, X21 has a non-negligible contribution to 

Xcff- 

Based on the parameters obtained from the fitting, we 
calculate the matrix elements function of the 

flux bias [Fig. 4(a)]. Remarkably, \gi2\ is more than two 
times larger than \goi\ at / = 0.5. In Fig. 4(b), we show 
the flux dependence of Xij' s i R Eq. (5). Only the three 
main components (xoij X12, X21) and XcS with N = 4 
are plotted. Reflecting the large \gi2\ and the fact that 
\ui2 - U) T \ ~ |w i - w r |, IX12I is much larger than |xoi| 
which corresponds to the dispersive shift estimated under 
the two-level approximation. Moreover, because we are 
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FIG. 4. (Color online) Flux bias dependence of \gtj\, Xij> 
and 2|x e fj|. (a) Calculated matrix elements of the cou- 
pling Hamiltonian H c between the two states \i) and \j), i.e. 
\gij\ — |{i|iT c |j)|. (b) Calculated \ij and Xeff- The three 
largest components are plotted here. The Xcff includes the 
terms up to j = 4 in Eq. ((5|. (c) 2|x ff| obtained from the 
measurement (blue dots), estimated from the energy band cal- 
culation (dashed red line), and calculated from Eq. (0) (solid 
green line). Experimental Xefi was obtained by Lorentzian fit 
to the measured spectra, and has errors of a few MHz. 



in the straddling regime, cJoi < w r < LJ12, all the three 
main components of x%j 's sum up in a constructive way in 
Eq. (5), which leads to the large | Xcff | ■ At /=0.5, 2|x e ff| 
is calculated to be 27rx72.7 MHz, which is consistent with 
the band calculation and the experiment. 

To further confirm the validity of our theoretical 
model, we measured the dispersive shift as a function of 
the flux bias, which is shown in Fig. 4(c). In this measure- 
ment, we obtained \XeB I m & different way from Fig. 3(a). 
We performed qubit spectroscopy as shown in Fig. 2(b), 
but with higher power of the qubit-control microwave. 
Due to the ac Stark shift, we observe side peaks around 
the main resonant peak of the qubit. Those peaks corre- 
spond to different photon numbers in the resonator Q, 
and their separation is equal to 2 |XefF | ■ As shown in the 
figure, the agreement between the experiment and the 
theory is fairly good. 

By exploiting the large dispersive shift, we measured 
the Rabi oscillations using the pulse sequence shown in 
the inset of Fig. 3(b). After a control pulse with the 
length of At is applied, we wait for i<j before we start 
acquiring the readout signal to calculate T. As shown 
in Fig. 3(b), we attained the contrast more than 90% 
after we experimentally optimized id to be 75 ns. We 



also estimate from the response time for the volt- 
age/current in the resonator 2k _1 = 2Q/u! r and the en- 
ergy relaxation time T\ of the qubit. We measured T\ 
at /=0.5 to be 680 ns (data not shown). The prod- 
uct (1 — exp(— 2t/ k)) ■ exp(— t/T\) shows its maximum 
at t =70 ns, which agrees well with the experimentally 
optimized t<\- Because 1 — exp(— t&/T\) amounts to nearly 
0.1, we attribute the 10% loss of the contrast to the en- 
ergy relaxation of the qubit during t&. 

In conclusion, we demonstrated strong coupling be- 
tween the flux qubit and the CPW resonator via a ca- 
pacitor. We observed a large dispersive shift of the cav- 
ity resonance and quantitatively accounted for its mag- 
nitude. The keys for the observed large dispersive shift 
are the large matrix element |<?i2| and the condition 
cjoi < Lo r < u\2 to realize the straddling regime. In ad- 
dition, we have demonstrated the Rabi oscillations with 
more than 90% contrast. Because the capacitive coupling 
for the flux qubit is advantageous in terms of the read- 
out backaction, the present result is an important step 
toward the high-fidelity single-shot readout of the flux 
qubit usin g J osephson bifurcation [3], [l4j or parametric 
amplifiers [151 ]. 
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bernetics" , the JSPS through its FIRST Program, and 
the NICT Commissioned Research. 



Appendix A: Derivation of the full Hamiltonian 

Figure 5 shows the equivalent circuit diagram of the 
device we measured. We assume that the larger Joseph- 
son junctions have capacitance Cj each, and the smaller 
one has aCj. In the diagram, L r , C r , and C c repre- 
sent the equivalent inductance and capacitance of the 
resonator, and the coupling capacitance between the res- 
onator and the flux qubit, respectively. The phase differ- 
ence across the larger Josephson junctions and the induc- 
tance L T are denoted as 8i (i=X, 2) and S r , respectively. 
The Lagrangian of the system is given by 




FIG. 5. Circuit diagram of the flux qubit capacitively coupled 
to the LC resonator. 
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where i?j = I <i>o/27r and / = $ ox /<I>o. This can be 
rewritten as 
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From this Lagrangian, we calculate the Hamiltonian 
of the system. We define the generalized momenta (p$ = 



acM), 

Pi \ / m'i -m' 3 m c \ / o\ \ 

P2 = -mj m' 2 -m c 5 2 , (A8) 

Pr / \ m c -m c / \8r / 



p = M8. (A9) 

By the Legendre transformation, "H = J^Pi^i ~ A we 
obtain the Hamiltonian 

H = ^p t M~ 1 p-U , (A10) 

where 

[/ = £j cos^i+^j cosS 2 +aEj cos(S 1 -5 2 +2nf) , (All) 



M = 



detM 
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m' 3 m[ — mr c m' x m' x — m 2 (m' 1 — m' 3 )m c 
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and 

det(M) = (m[m 2 -m 3 2 )m' r +(2m 3 -m[-m' 2 )m 2 c . (A13) 

The Hamiltonian consists of three parts, the qubit part 
Hq, the resonator part "H r , and the coupling part W c . 

"H = + "Hr + He- (A14) 
The qubit part is given by 
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and By using 
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P = C c /Cj, (A18) "~ 

7 = C c /C r . (A19) "H q is written as 
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The resonator part is given by 
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From this, we get the loaded resonant frequency u> T , 
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Using creation and annihilation operators, we write the 
resonator part of the Hamiltonian as 



"H r = h(jj t {a^a + — ) . 
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Finally, the coupling part is given by 
m' 



where 
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Appendix B: Method 

1. Experimental setup 

Figure 6 shows an experimental setup diagram for our 
measurements. Qubit control and readout microwave 
pulses are generated by mixing the continuous microwave 
with pulses of 50 MHz IF frequency generated by DACs 
developed by Martinis group at UCSB [16|. The pulses 



are sent through the input microwave lines each of which 
include attenuators of 42 dB in total. For the read- 
out line, the microwave pulses are further attenuated by 
20 dB, and are routed to the resonator through a circu- 
lator to separate the input and output (reflected) waves. 

The output line for readout includes a band-pass filter 
(9-11 GHz), a low-pass filter (/ c =12.4 GHz), and three 
isolators with ~20 dB isolation each (9-11 GHz) in order 
to prevent noise/blackbody radiation coming from the 
cryogenic HEMT amplifier. The reflected signal is am- 
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plified by the amplifier and a room-temperature amplifier 
with a total gain of ~66 dB, and mixed with a local os- 
cillator at an I/Q mixer down to the IF frequency. The I 
and Q quadratures are sampled at I GS/s by a digitizer. 
The data is averaged 6.5 xlO 4 times for each measure- 
ment point. 

For the frequency-domain measurement (Fig. 2), we 
replace the readout and the control circuits with a vector 
network analyzer and a continuous microwave generator, 
respectively. 

The excitation energies of the flux qubit are controlled 
by the local flux bias $ which is generated by a supercon- 
ducting coil wound around the sample package mounted 
inside the cryogenic magnetic shield. 



2. Spectroscopy data fitting 

Figure 7(a) and (b) show the fitting results of the mea- 
sured spectra (the same as Fig. 2(a) and (b)) with the 
energy levels obtained from the Hamiltonian (Eq. IA14[) . 
The spectra are well fitted using 1 1 charge states for each 
islands of the flux qubit and the number of levels for the 
resonator is truncated by 5. The following circuit pa- 
rameters are extracted by the fitting of the energy lev- 
els: E 3 /h=U8A GHz, E c /h=3.268 GHz, a=0.6106, and 
C c =4.079 fF. For comparison, the same data as in (a) and 
(b) are shown in (c) and (d), respectively. 



3. Normalization of measured signals 

Figure 8(a) shows the reflection coefficient V measured 
with the pulsed microwave and the heterodyne setup 
shown in Fig. 6. The data is represented in the po- 
lar coordinates for the frequency range from 10.50 to 
10.75 GHz. The point A (B) indicated by the green (pur- 
ple) arrow gives T at ui p /2ir = uj r /2ir =10.656 GHz when 
the qubit is in |0) (|1)) state. When we measure the Rabi 
oscillations, namely, fix lo p at lo v and vary the length of 
the qubit control pulse, observed T goes back and forth 
between the point A and B. The dashed line in Fig. 8(a) 
represents a locus of the signal of the Rabi oscillations. 
Thus, the maximum amplitude of the Rabi oscillations 
is given by |AB|. The point C indicated by the light 
blue arrow is T at an off-resonance of w p /27r=10.52 GHz, 
which is the reference data used to normalize the mea- 
sured data. We define the normalized reflection coeffi- 
cient T' = (r -^r on )/|r o gf - r on |, where r on and r o ff 

correspond to A and C, respectively, and |r o g — r on | 
corresponds to the expected maximum amplitude of the 
Rabi oscillations. Figure 8(b) depicts V represented in 
the polar coordinates. The points A, B, and C in (a) are 
displaced to A' (the origin), B', and C, respectively. |r'| 
is equal to zero at uj v = uj r when the qubit is in |0) state, 

and is equal to |A'B'| corresponding to the normalized 
maximum amplitude of the Rabi oscillations ^0.9 when 
a 7r-pulse is applied. 
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FIG. 6. Experimental setup diagram for the heterodyne detection using the pulsed microwave. 




FIG. 7. Spectra of the coupled system. The dashed curves in (a) and (b) represent the fit to the energy levels calculated from 
the Hamiltonian, Eq. (|A14[) . For comparison, the same data as in (a) and (b) are shown in (c) and (d), respectively. The data 
in (a) [(c)] is embedded in (b) [(d)] in the blue box. There is no data in the regions filled with the meshes. 
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FIG. 8. Normalization of the reflection coefficient T. (a) T represented in the polar coordinates for the frequency range from 
10.50 to 10.75 GHz. The point A (B) corresponds to F at u p /2tt = ^ r /27r=10.656 GHz when the qubit is in |0) (|1)) state. 
The point C is F at an off-resonance of u> p /2n =10.52 GHz, which is a reference data used to normalize the measured data . 
(b) Normalized reflection coefficient T' represented in the polar coordinates. The points A , B and C in (a) are displaced to A' 

(the origin), B', and C', respectively. The normalized maximum amplitude of the Rabi oscillation corresponds to |A'B'| and it 
is ~0.9. 



